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How does Fermi surface develop in cuprates upon doping of a parent Mott insulator, does it 
consist of large barrels or small pockets, which of them is responsible for superconductivity and 
what is a role of the pseudogap? Those are actively debated questions, important for understanding 
of high temperature superconductivity. Here we analyze doping dependence of interlayer tunneling 
in cuprates. We observe that with decreasing doping the supercurrent is rapidly decreasing, but the 
quasiparticle resistance at high bias remains almost unchanged. This indicates that Cooper pair and 
quasiparticle currents originate from different parts of Brillouin zone: Cooper pairs are residing only 
on small pockets, which are progressively shrinking with decreasing doping, but the quasiparticle 
current is integrated over the full length of barrels, which are only weakly doping dependent. The 
expanding pseudogap areas along the barrels do not contribute to pair current. This provides direct 
evidence for nonsuperconducting origin of the pseudogap. 


Fermi surface in metals occurs at an intersection of a 
conduction band with a chemical potential. Insulators 
do not have Fermi surface because the chemical poten¬ 
tial lies in the band gap region. High temperature super¬ 
conductivity in cuprates appears upon doping of a Mott 
insulator. One of the key questions is how Fermi surface 
develops with doping m- According to Luttinger theo¬ 
rem Fermi surface area should be proportional to doping 
p. However, for cuprates the photoemission edge, seen by 
angular-resolved photoemission spectroscopy (ARPES), 
does not grow gradually with doping but forms large bar¬ 
rels with an area oc (1 +p) already at the lowest doping 
[H]. However, strictly speaking only nodal parts of the 
barrels (Fermi arcs) are representing true Fermi surface 
because anti-nodal parts are gapped HIS]. The arcs size 
is growing linearly with doping [MSI- Luttinger theo¬ 
rem can be satisfied assuming that arcs represent parts 
of small Fermi pockets gmniiii]- Although existence 
of small pockets was confirmed by quantum oscillation 
experiments [T1HI3, their position in Brillouin zone [T], 
connection to barrels and significance for high-T^ super¬ 
conductivity remains unclear m- 

The role of a normal state pseudogap (PG) is another 
related issue [13 [18]. Similarities between the PG and 
the superconducting gap (SG) have led to an assumption 
of precursor superconductivity origin of the PG |19j . In 
this case anti-nodal PG parts of barrels should contain a 
major part of the superconducting condensate. However, 
there are also arguments in favor of competition of the 
two co-existing gaps 0 m [13 HOMS]- Discrimination 
between SG and PG is particularly difficult at T < Tc, 
when the whole barrel is gapped. So far it is not possi¬ 
ble to conclude whether superconductivity is originating 
from large barrels or small Fermi pockets. 

Here we analyze doping dependence of interlayer tun¬ 
neling characteristics of small Bi 2 Sr 2 Cai_a;Ya;Gu 208 -i -5 
[Bi(Y)-2212] intrinsic Josephson junctions. We utilize the 
ability of superconducting tunnel junctions to indepen¬ 
dently probe Gooper pair and single quasiparticle (QP) 
currents. We observe that the high-bias QP resistance re¬ 


mains almost doping independent, implying that QP cur¬ 
rent is originating from weakly doping dependent barrels. 
To the contrary, the supercurrent is rapidly decreasing 
with decreasing doping, indicating that Cooper pairing 
occurs only on small Fermi pockets that progressively 
shrink with decreasing doping. The antinodal parts of 
barrels, which grow with underdoping, do not contribute 
to supercurrent. This directly proves that the PG is not 
due to Cooper pairing. We present numerical calcula¬ 
tions that support our conclusions. 

Intrinsic Josephson junctions are naturally formed in 
Bi-2212 single crystals m- Atomic scale of such junc¬ 
tions leads to a large capacitance and a quality factor Q ^ 
10^ [28]. Current-Voltage {I-V) characteristics of junc¬ 
tions with Q 3> I acquire a hysteresis with zero-voltage 
and resistive branches corresponding to Cooper pair and 
QP tunneling, respectively, which allows 

independent analysis of pair and QP transport. 

We study Bi(Y)-2212 single crystals from the same 
batch with an optimal Tc{OP) ~ 95 K. Doping level p 
was changed by annealing at T = 600° C and was esti¬ 
mated using an empirical expression Ti,{p) = Tc{OP)[l — 
82.6(p - 0.16)2] [Hj. It corresponds to an onset of su¬ 
perconductivity at p = 0.05 at the insulator-to-metal 
transition and to the optimal doping (OP) at p = 0.16 
(holes per Cu). Micron-size mesa structures, contain¬ 
ing A^ = 7 — 12 junctions were made using micro/nano¬ 
fabrication techniques. Details of sample fabrication and 
characterization can be found in Refs. E5] . 

Figure [^ summarizes temperature and doping depen¬ 
dencies of interlayer tunneling. Fig. (a) shows T- 
dependencies of I-V's for a moderately underdoped (UD) 
mesa p ~ 0.14. At T < Tc a sum-gap kink appears at 
V = 2/1^1 e (per junction) followed by Ohmic, and almost 
T-independent tunnel resistance i?„. The constancy of 
Rn is a fundamental consequence of electronic state con¬ 
servation [^. Spectroscopic features are better analyzed 
using dI/dV{V) characteristics, shown in Fig. 0(b). It 
is seen that the sum-gap peak shifts to lower voltages 
and decreases in amplitude with increasing T and van- 
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FIG. 1. (Color online) Temperatnre and doping dependence of intrinsic tunneling characteristics, (a) I-V and (b) dl/dViV) 
characteristics of a small moderately nnderdoped mesa at different T. Inset demonstrates crossing of dl/dViV) curves at 
T > Tc. (c) Temperature dependencies of the superconducting peak and the pseudogap hump voltages for the same mesa, (d) 
Comparison of normalized intrinsic spectra for slightly and strongly underdopoed mesas, (e) and (f) Comparison of normalized 
I-V curves at different doping (e) in the superconducting and (f) in the normal states. 


ishes at Tc- At T > Tc the zero-bias conductance remains 
suppressed as a consequence of persisting PG. The corre¬ 
sponding missing states are expelled into a broad hump 
at higher voltages. With increasing T the zero bias mini¬ 
mum fills-in and the hump is deflated in a state conserv¬ 
ing manner ED 12 a, so that the dl/dV{V) characteristics 
at different T cross at one point, as shown in the inset 
ED US- Fig. a (c) shows T-dependencies of the super¬ 
conducting peak and the pseudogap hump. It is seen that 
the SG decreases rapidly upon approaching the in a 
BGS (Bardeen-Cooper-Schrieffer) manner ED- The PG 
persists in a broad T range both below and above Tc- 

Fig. 0 (d) shows differential conductivity for a slightly 
p = 0.143 and strongly p = 0.123 UD mesas at low T. It 
is seen that for a strongly UD mesa p = 0.123 the hump 
coexists with the peak at T ^ Tc EOl ED , the peak is 
strongly suppressed, the hump is enhanced and has even 
larger amplitude than the peak. With increasing doping 
the peak height is increasing and the hump is decreasing 
both in height and voltage. With further increase of dop¬ 
ing the hump is buried under the peak, as seen from Fig. 

(d) for p = 0.143. The hump, however, is uncovered at 
elevated T as the peak shifts to lower voltages ED- 

Fig. 0 (e) and (f) represent current density J vs. volt¬ 


age per junction characteristics for different doping levels 

(e) at low T and (f) above Tc. The following main fea¬ 
tures, which will be in focus of our discussion, are seen: 

(i) The critical current. Multiple branches at low bias 
appear due to one-by-one switching of junctions from the 
superconducting to the resistive state EqIEt]. The am¬ 
plitude of the branches represent the Josephson critical 
current density Jc, marked in Fig. [^(e). Jc rapidly de¬ 
creases with decreasing doping. As a result, branches, 
which are very pronounced for the overdoped (OD) mesa 
p = 0.182, are hardly visible (on this scale) for the most 
UD p = 0.113 mesa. Reduction of Jc reflects reduction of 
the superconducting condensate with decreasing doping. 

(ii) The sum-gap kink amplitude. AJqp represents a 
number of QP states within the superconducting gap, 
which are subjected to pairing. Therefore, AJqp and 
Jc should be directly connected. Indeed, from Fig. 

(e) it is seen that AJqp is also rapidly decreasing with 
decreasing doping so that the kink becomes poorly visible 
for the most UD mesa. This indicates a reduction of the 
number of electronic states subjected to Gooper pairing. 

(iii) The pseudogap is most clearly visible in the normal 
state when superconducting features are gone. From Fig. 

(f) it is seen that the I-V of a strongly UD mesa p = 
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FIG. 2. (Color online) Numerical modelling of intrinsic tunneling characteristics, (a) Sketch of a quarter of the Brillouin zone. 
(b)-(d) Angular distribution along the barrel from an antinodal to a nodal point of (b) the gaps, (c) the quasiparticle damping 
factor and (c) the tunneling transmission probability, (e)-(g) Evolution of calculated dI/dV{V) curves upon varying of (e) the 
pocket size, (f) The pseudogap energy and (g) the superconducting gap. Arrows indicate positions of peaks and humps. 


0.113 is nonlinear at T > Tc. The current is suppressed 
below a threshold voltage AVpc- AVpc decreases with 
increasing doping and disappears for slightly OD mesas, 
as seen from an almost Ohmic I-V at p = 0.182. 

(iv) The tunnel resistance. From Fig. (e) it is seen 
that the high-bias Rn is initially increasing with decreas¬ 
ing doping, but for UD mesas the differential resistance 
becomes almost doping independent. Since tunnel re¬ 
sistance depends on the QP DoS it is expected that 
should increase with decreasing carrier concentration and 
saturation at low doping requires explanation. 

To gain a better understanding of experimental char¬ 
acteristics we performed numerical analysis taking into 
acconnt a symmetry of the order parameter and a topol¬ 
ogy of the Fermi surface, as shown in Fig. (a). Follow¬ 
ing recent studies [lOHI^ , we assume that central parts of 
large barrels represent small Fermi pockets with an angu¬ 
lar size if Arc- Fig- [^(b) represents angular dependencies 
of the gaps along the barrel. In the normal state only 
antinodal PG parts of barrels are gapped. At T < Tc a 
d-wave SG is opening. We assume that the antinodal re¬ 
gions acquire a combined gap Acomb = \/^pg + ^SG- 
The exact scenario in the antinodal region is not very im¬ 
portant because the main difference between barrel and 
arc regions is in the QP damping factor F, which is small 
at the pocket and large at the barrel, as shown in Fig. 
(c). This makes QP’s ill-defined in the PG region [2]. Fi¬ 
nally, we have to take into account angular dependence of 
the transmission probability Tr{ip). For non-directional 
tunneling Tr{ip) = const. However, c-axis transport in 


cuprates should be dominated by antinodal regions. The 
corresponding directional Tr{ip) is shown by the dashed- 
dotted line in Fig. [^(d). We have found that the experi¬ 
mental data is best fitted using a semi-directional Tr((p), 
which is finite at nodal regions, as shown by the solid line 
in Fig. [^(d). Numerical simulations presented in Figs. 

(e-g) are made for the semi-directional case. Details of 
calculations can be found in the Supplementary [51] . 

Fig. (e) demonstrates variation of dI/dV{V) upon 
changing of the arc size. The case pArc = 90° corre¬ 
sponds to the absence of the PG. In this case the spec¬ 
trum contains a single superconducting peak at eV = 
2A5(3(0°) with the height and the shape similar to that 
for OP mesas. As the arc shrinks, the amplitude of the 
peak is rapidly decreasing and the PG hump is growing 
at eV = 2Acomh(0°). This is similar to evolution of ex¬ 
perimental curves with decreasing doping in Fig. [I](d). 

Fig. (f) demonstrates variation of dl/dViy) upon 
changing the PG energy. When Apc is significantly 
larger than Asg both the peak and the hump are well 
defined. However, as Apq starts to decrease the hump is 
moving towards the peak and is eventually buried under 
the peak. This is similar to evolution of experimental 
curves with increasing doping, see Figs, [^(b) and (d). 

Fig. (g) shows variation of dljdV (P) upon changing 
of AsGt mimicking T-variation shown in Fig. [^(b). Note 
that the hump moves to lower voltages with decreasing 
Asg because in our case it occurs at the combined gap. 
We observe a similar shift of the hump at T < Tc in 
experiment, as shown in Fig. (c). We conclude that 


















4 




FIG. 3. (Color online) Doping dependence of (a) the pseudogap, (b) high-bias resistivities, (c) the critical current density and 
the amplitude of the sum-gap kink AJqp. (d) Calculated dependence of the critical current from Fermi pockets, as a function 
of ifiArc for non-directional, semi-directional and directional tunneling (dashed, solid and dashed-dotted lines), (e) Comparison 
of scaled Ic and AlQpj2. They represent the amounts of Cooper pairs and QP states subjected to pairing, respectively. Lines 
(top-right axes) represent normalized curves from panel (d). The reduction of 7c and AIqp with decreasing doping is consistent 
with proportional shrinkage of the Fermi pockets containing the superconducting condensate, (f) The solid line represents a 
deduced doping dependence of the pocket size. Symbols represent corresponding data obtained by other techniques. 


there is a good overall agreement between the experimen¬ 
tal data and the considered model. Therefore, we employ 
this model for a quantitative analysis of data. 

Figure represents a summary of doping dependence 
of intrinsic tunneling characteristics. Fig. |^(a) shows the 
PG energy obtained from half of the PG hump voltage. 
It is linearly decreasing with increasing doping and tends 
to vanish at the quantum critical point slightly above 
Pc = 0.19, in agreement with Refs. [siin]. A topological 
barrel-pocket transition should occur at this point |12) . 

In Fig. (b) we analyze doping dependence of high 
bias resistivity. Solid and open symbols represent dc and 
differential (ac) values, respectively. From Fig. 0 (f) it 
is seen that for OD mesas the I-V is Ohmic and the 
two resistances coincide. The /-P’s of UD mesas become 
non-linear due to appearance of the PG and the two resis¬ 
tances become different. The linear doping dependence 
of the dc resistivity reflects the corresponding behavior of 
the PG, Fig. 0(a). However, the ac resistivity becomes 
almost doping independent at low doping, as seen from 
Fig. 0 (e). Simulations provide a clarification of different 
behavior of dc and ac resistivities. From Figs. 0 (e) and 


(f) it is seen that the high-bias dl/dV is independent of 
If Arc and Apc- This occurs because the PG does not 
change the total amount of states but just redistributes 
them. At a sufficiently high voltage eV ^ Apc all QP 
states along the barrels, including the gapped anti-nodal 
parts, contribute to the QP current. Therefore, the weak 
doping dependence of ac-resistivity is a consequence of 
the weak doping dependence of barrels laa. 

Fig. |3 (c) shows doping dependence of Jc and AJqp It 
is seen that unlike i?„ both Jc and AJqp decay rapidly 
with decreasing doping. Behavior of Jc{p) is in agree¬ 
ment with previous reports O A qualitative 

difference of doping dependencies of the Gooper pair and 
the high-bias QP transport is the main new observation 
of this work. As explained above, the weak doping depen¬ 
dence of Rn suggests that the QP current is integrated 
over the full length of barrels, which are weakly doping 
dependent. Therefore, a rapid decrease of the Gooper 
pair current with decreasing doping indicates that the su¬ 
perconducting condensate does not reside along the full 
length of barrels, but occupy a progressively smaller frac¬ 
tion, as expected for small Fermi pockets. 
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To understand how I,, depends on the pocket size we 
performed corresponding calculations, shown in Fig. 

(d) (see Ref. [34]). In the non-directional case, 

Ic = [1 - COs((/?Arc)] ■ (1) 

In the absence of the PG, tpArc = 90°, IcRn = Ao/e. 
With decreasing doping Fermi arcs shrink and Ic de¬ 
creases. The Ic vanishes when arcs collapse (pArc 0. 
The supercurrent and the high-bias QP resistance have 
different doping dependencies because supercurrent is 
measured at zero voltage. Consequently, pair tunneling 
occurs only between true (ungapped) Fermi surfaces, i.e., 
Fermi pockets, which shrink with decreasing doping. To 
the contrary, QP tunneling occurs at finite bias and at 
eV A PC it collects all QP states along the barrels. 

From Fig. |^(c) it is seen that AIqp and Ic are chang¬ 
ing in a correlated manner. For a d-wave tunneling in the 
absence of the PG the kink amplitude should be twice the 
critical current AIqp = 2Ao/ei?„ = 2Ic. In Fig. [^(e) we 
show the correspondingly scaled quantities as a function 
of doping. They were obtained using the high bias dif¬ 
ferential resistance as i?„ (open and solid symbols were 
obtained using different criteria for estimation of i?„ and 
A). It is seen that Ic and AIqp/2 merge together, con¬ 
firming that both are indeed related: Ic is measuring the 
amount of Cooper pairs and AIqp the amount of QP 
states within the gap, subjected to pairing. Importantly, 
the scaling indicates that both quantities originate from 
the same parts of the Brillouin zone ( i.e., Fermi pock¬ 
ets). The reduction of both Ic and AIqp with decreasing 
doping reflects the shrinkage of the pocket size. 

In order to estimate variation of <pArc with doping, in 
Fig. i (e) we compare experimental data with theoretical 
curves from Fig. [^(d). Here we assumed that ip Arc = 90° 
at the onset of the PG, p = 0.192, and that ipArc = 
0° at the insulator-to-metal transition, p = 0.05. The 
agreement is remarkable taking into account that there 
is no fitting other than adjustment of the vertical scale. 

Figure (f) summarizes our main result. A solid 
line represents the deduced angular size of the Fermi 
pocket/arc as a function of doping. It shrinks linearly 
with decreasing doping, consistent with previous reports 
by other techniques liiHiiniiii], shown for comparison in 
the figure. An important new aspect of our work is that 
we could discriminate between Cooper pair and quasi¬ 
particle transport. We conclude that the QP current at 
high bias is originating from the full area of large bar¬ 
rels, leading to weak doping dependence of i?„. However, 
the Cooper pair current originates only from small Fermi 
pockets, which shrink with decreasing doping, leading 
to a rapid decrease of both the critical current and the 
sum-gap kink amplitude. The pseudogap parts of barrels, 
which grow bigger with underdoping, apparently do not 
contribute to supercurrent. This is a direct evidence for 
non-superconducting origin of the pseudogap in cuprates. 


Technical support from the Core Facility in Nanotech¬ 
nology at SU is gratefully acknowledged. 


SUPPLEMENTARY INFORMATION 

The supplementary material provides additional infor¬ 
mation about numerical calculations of interlayer tun¬ 
neling characteristics and analytical results for variation 
of IcRn as a function of the arc size for non-directional, 
directional and semi-directional energy and momentum 
conserving tunneling between d-wave superconductors. 
We assumed that the supercurrent is originating only 
from nodal arcs and single quasiparticle current from the 
full length of large barrels. 


Calculation of interlayer tunneling characteristics 


Proper calculation of interlayer tunneling character¬ 
istics requires accurate integration of tunneling current 
over all initial and final states, taking into account the 
band structure, angular dependence of the energy gap 
A(y)) and the transmission probability Tr(ipi,ip 2 ) for 
tunneling between the initial and final states with the 
momentum angles ipi^ 2 y where ip is the angle between the 
momentum of the QP and the principle axis of the Bril¬ 
louin zone. We consider the case of elastic tunneling, in 
which the energy of electrons is conserved. In this case 
the QP tunneling current between two adjacent super¬ 
conducting layers can be written as: 


r2Tr 


I = A 


dpi 

27T 




( 2 ) 


Tr\puI’ 2 )N{E, pi)N{E + eU, P 2 ) [f{E) - f{E + eV )], 


where E is the energy of the QP with respect to the 
chemical potential, N{E, pi) and N{E -\- eU, P 2 ) are the 
corresponding QP DoS in the initial and final states: 


N{E,ip) = Y(0)3fi 


E - ir((p) 


(3) 


Here r(<p) is the angular dependent quasiparticle damp¬ 
ing factor (the inverse QP lifetime). In the absence of 
the PG, analysis of various tunneling scenarios for coher¬ 
ent ipi = p 2 and incoherent pi ^ ip 2 tunneling between 
d-wave superconductors can be found in the supplemen¬ 
tary material to Ref. [29] and in Ref. [35]. Here we 
will focus on analysis of the pseudogap effect for the case 
of coherent ipi = p 2 and elastic (momentum and energy 
conserving) tunneling for a model of a “remnant” Fermi 
barrel with a PG in antinodal and a Fermi-arc in the 
nodal parts, as shown in Fig. 2 (a). We assume that 
Cooper pair tunneling occurs only in the Fermi-arc re¬ 
gions, but the QP tunneling occurs along the whole bar¬ 
rel. All the characteristics are presented for a quarter of 
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Brulloin zone and are symmetrically reflected for other 
three quarters of the zone. 

Fig. 2 (b) represents angular dependence of the gaps. 
In the normal state only PG is present in the antinodal 
parts of the barrels in the angular interval from 0 to 
^PG = 7r/4 — ipArcl‘2.- We assumed the following angular 
dependence, 

^pg{^) = Apg(O) cos > (4) 

shown by the dashed-dotted line in Fig. 2 (b). 

At T < Tc there is a superconducting gap with a d- 
wave symmetry, 

Asg(¥?) = Aocos(2(^), (5) 

shown by the dashed line in Fig. 2 (b). We assume that 
in the PG region the two gaps form a combined gap 

Acombi^p) = ^A^q(v 3) -P A|g(v5), (6) 

shown by the solid line in Fig. 2 (b). This assumption is, 
however, not critical because antinodal parts of the barrel 
do not contribute to a sharp sum-gap peak in dl/dV 
because of a large QP damping in the PG region [5], 
as shown in Fig. 2 (c). We assumed that within the 
PG region the F pq is varying in a similar manner to 
Eq. Q with the maximum rpG(O) ~ lOmeV. In the 
arc region the QP damping is two orders of magnitude 
smaller Psg = O.ImeV. Due to the large difference in 
F sharp sum-gap kink/peak features in I-V and dl/dV 
originate solely from the arc region and correspond to 
the superconducting gap. Eq. ([^ simply represents a 
comfortable way to connect the PG and SG regions by a 
continuous line. 

For the angular dependence of the transmission prob¬ 
ability Tr{(p) we considered three scenarios: 

(i) Non-directional tunneling Tr((p) =const, as shown 
by the dashed line in Fig. 2 (d). 

(ii) Directional tunneling with maximum in antinodal 
and zero transmission in nodal regions IMti] . Tr{ip) oc 
[cos(A:a;) — cos(A; 2 ^)]^ oc [cos(7rsin((^)) — cos(7rcos((^))]^, as 
shown by the dashed-dotted line in Fig. 2 (d). This ex¬ 
pression is well approximated by the function cos(2(p)^, 
which we will use for analytical calculations below. 

(hi) Semi-directional tunneling, which is an average 
of non-directional and directional cases, Tr{ip) cx 1 -I- 
[cos(7rsin(i^)) — cos(7r cos((^))]^/4, as shown by the solid 
line in Fig. 2 (d). 

Finally, the QP density of states N{ip) is likely to vary 
along the barrel. However, in the absence of a confident 
knowledge of the angular dependence of DoS we just as¬ 
sumed it to be constant N{(p) = N{Q). The angular 
variation of DoS would have the same affect as angular 
variation of the transmission coefficient. 


Figure shows calculated I-V characteristics in case 
of semi-directional tunneling, corresponding to dl/dV 
curves in Fig. 2 (e-g). Figs. (a-c) demonstrates vari¬ 
ation of I-V's upon changing (a) the arc size, (b) the 
PG energy and (c) Asgi while keeping other parameters 
constant. Note that the high bias resistance remains the 
same, irrespective of the PG. This occurs because in our 
case the PG is a state-conserving gap, see Eq. There¬ 
fore the same current is recovered upon integration over 
the full barrel at high enough bias, irrespective whether 
there is a PG or not. Note that curves in Figs. |^(b) and 
(c) calculated for a large F pq (0) = 20 meV resemble the 
experimental characteristics for UD mesas with a finite 
offset voltage AVpc, see Figs. 1 (e) and (d). 

Dependence of the critical current on the arc size 

For a junction made of s-wave superconductors IcRn is 
independent of the transmission coefficient and is given 
by Ambegaokar-Baratoff expression, 

ttA 

hRn{T^T,) = —. (7) 

2e 

For coherent tunneling in a junction made of d-wave 
superconductors, Eq. Q is still valid for a specific an¬ 
gle [37],i.e. for 6Ic{p)6Rn{<p), where Sldp) and 6Rn{p) 
are contributions to the critical current and resistance 
from Gooper pairs and QP’s with momentum in the di¬ 
rection tp. Thus, SIc{<p) = TTA{(p)/2eSRn{>p)- From Eq. 
([^ l/dRn{p) = AeN{0)'^Tr{ip)^ Sipl2n. Introducing a 
quantity Ri = AeN{0)'^ for resistance at unit transmis¬ 
sion Tr = 1, we can write expressions for the total critical 
current and the normal resistance, taking into account d- 
wave angular dependence of the gap A{(p) = Aq cos(2(p): 

4 

IcRi = — / Tr{ipf‘cos{2ip)dip (8) 

JIT/‘i-tpAral'i 

R,R-^ [ ' Tr^pfdp (9) 

^ Jo 

Here we assumed that Cooper pair current is originating 
only from the arc and single QP current is accumulated 
over the full length of the barrel. 

After integration of Eqs. Q and we obtain : 

IcRn = ^ [1 - COs(:/JArc)] , (10) 

for the nondirectional case Tr{ip) oc I. For tpArc = '^/‘I 
it reduces to a known expression IcRn — Ag /em- 


COs((/^^j’c) “1“ Arc) _ COs((p^7-c) 5 

15 3 5 
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FIG. 4. Numerical modelling of intrinsic tunneling characteristics. Evolution of calculated I-V curves upon varying 
of (a) the pocket size, (b) The pseudogap energy and (c) the superconducting gap. Note that the high-bias resistance remains 
unchanged because it is integrated over the same barrel. 


for the directional case, Tr{(p) oc cos(2(/3)^, and 


— - 4 COs{(pArc) + o “ 7 C0s((/:7Arc)^ , 

15 3 5 

for the semi-directional case, Tr{(f) oc 1 -b cos(2(/?)^. The 
corresponding curves IcRni^PArc) are shown in Fig. 3 
(d). In all cases Ic vanishes as arcs collapse (pArc —t 0 
because we assumed that the supercurrent is originating 
only from arc regions. The i?„ is independent of (pArc 
because the QP current is integrated over the whole bar¬ 
rel, which we assumed to be unchanged. As described 
in the manuscript, such model provides a qualitative ex¬ 
planation of different doping dependencies of the critical 
current density and the high-bias differential resistance. 
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